题目 城市 表层 土壤 重金 属 污染 分 析 


作 ж: BAR 
带 队 老师 : EEE 
郑州 市 57 中 


摘要 : 随 着 人 类 活动 对 城市 环境 质量 日 益 突 出 的 影响 ， 对 城市 土壤 地 质 环境 异 常 的 查证 


逐渐 成 为 人 们 关注 的 焦点 。 本 文 利用 采样 点 的 数据 建立 了 数学 模型 对 城市 环境 的 质量 做 
出 评价 与 判断 。 针 对 题目 要 求 ， 首 先 绘制 出 重金 属 分 布 图 与 海拔 等 值 线 图 ， 给 出 重金 属 
在 该 城区 的 空间 分 布 。 然 后 以 《土壤 环境 质量 标准 》(GB15618-1995) 为 基本 标准 给 出 
污染 评价 。 通 过 图 的 对 比 和 污染 评价 得 出 重金 属 污 染 原 因 ， 传 播 特征 及 污染 源 位 置 。 最 
后 对 模型 作出 评价 。 

问题 一 : 采用 三 次 多 项 式 插值 的 方法 绘制 重金 属 分 布 图 和 海拔 等 值 线 图 ， 给 出 重金 
属 空间 分 布 。 然 后 采用 单 因 子 污染 指数 和 尼 梅 罗 污 染 综 合 指数 给 出 污染 评价 ， 得 出 生活 
区 、 交 通 区 、 工 业 区 为 重 污染 ， 公 园 绿地 区 为 中 度 污 染 ， 山 区 为 轻 度 污染 的 结果 ， 其 中 
Cu 和 Hg 的 污染 最 为 严重 。 

问题 二 : 根据 污染 物 迁 移 理论 结合 分 布 图 和 污染 的 分 区 评价 ， 参 考 重金 属 用 途 ， 分 
元 素 分 区 域 综 合 考 虑 ， 得 到 污染 的 原因 主要 为 工业 生产 、 生 活 垃 圾 、 燃 料 的 燃烧 、 汽 车 
尾气 等 。 

问题 三 ， 根据 每 种 重金 属 在 该 城区 空间 分 布 的 不 同 特征 ， 分 析出 Cr, Cu, Не. Ni 
呈 富 集 状 态 ，As、Cd 随 水 流 、 空 气 等 因素 传播 ，Pb、Zn 的 传播 受 人 为 因素 影响 较 大 。 
对 于 空间 分 布 广泛 的 重金 属 ， 建 立 赋 权 值 的 重心 选 址 模型 ， 采 用 穷 举 法 ， 得 出 重金 属 污 
染 源 的 大 概 位 置 ， 进 行 数据 筛选 ， 再 次 应 用 重心 选 址 模型 ， 得 到 污染 源 的 位 置 。 对 于 集 
中 分 布 的 重金 属 ， 直 接 进 行 数据 筛选 ， 通 过 重心 选 址 模型 确定 每 种 重金 属 元 素 的 污染 源 
位 置 坐 标 。 

问题 四 : 对 三 次 多 项 式 插 值 和 重心 选 址 模型 作出 评价 ， 并 分 别 运 用 克 里 金 插值 和 反 
演 模型 进行 改进 。 对 尼 梅 罗 污 染 综 合 指数 评价 法 赋 权 值 进行 改进 。 如 果 收 集 到 不 同 深度 
的 土壤 检测 数据 (重金 属 元 素 含量 、PH 值 。 土 壤 含水 率 等 )， 可 以 建立 推 流 迁 移 模 型 ， 
得 出 污染 时 间 的 变化 规律 以 及 污染 源 的 具体 位 置 。 

关键 词 ， 三 次 多 项 式 插值 ” 单 因 子 污染 指数 ” 尼 梅 罗 污 染 综合 指数 ”重心 选 址 


а. 


模型 


| Аы мен ана ААН 4 
РИ А 4 
ДЕС ЛЕН ИЛ Т Т ТТ 5 
Д: TH  еееееЗеЗЗее- 5 
Dl ra таа 6 
RE sd TREE 25 
Md aga 26 


城市 表层 土壤 重金 属 污染 分 析 
1 问题 重 述 

随 着 城市 经 济 的 快速 发 展 和 城市 人 口 的 不 断 增 加 ， 人 类 活动 对 城市 环 
境 质量 的 影响 日 显 突 出 。 对 城市 土壤 地 质 环境 异常 的 查证 ， 以 及 如 何 应 用 
查证 获得 的 海量 数据 资料 开展 城市 环境 质量 评价 ， 研 究 人 类 活动 影响 下 城 
市 地 质 环 境 的 演变 模式 ， 日 益 成 为 人 们 关注 的 焦点 。 

按照 功能 划分 ， 城 区 一 般 可 分 为 生活 区 、 工 业 区 、 山 区 、 主 干道 路 区 
及 公园 绿地 区 等 ， 分 别 记 为 1 类 区 、2 类 区 、……… 、5 类 区 ,不 同 的 区 域 环 
境 受 人 类 活动 影响 的 程度 不 同 。 

现 通过 对 某 城市 城区 土壤 地 质 环 境 进行 调查 ， 我 们 要 通过 数学 建 模 来 
完成 以 下 任务 : 

(1) 给 出 8 种 主要 重金 属 元 素 在 该 城区 的 空间 分 布 , 并 分 析 该 城区 内 不 
同 区 域 重 金属 的 污染 程度 。 

(2) 通过 分 析 数 据 ， 说 明 重 金属 污染 的 主要 原因 。 

(3) 分 析 重 金属 污染 物 的 传播 特征 , 由 此 建立 模型 , 确定 污染 源 的 位 置 。 

(4) 分 析 所 建立 模型 的 优 缺 点 ， 为 更 好 地 研究 城市 地 质 环境 的 演变 模 
式 ， 另 外 搜集 一 些 信 息 ， 简 述 建 立 模型 和 解决 问题 的 方法 。 

2 问题 分 析 

城区 按 功 能 划分 为 生活 区 、 工 业 区 、 山 区 、 主 干道 路 区 及 公园 绿地 区 ， 
不 同 区 域 环 境 受 人 类 活动 影响 的 程度 不 同 ， 根 据 题目 所 给 数据 分 别 建立 模 
型 求解 。 

问题 一 、 关 于 重金 属 的 空间 分 布 ， 题 目 所 给 信息 中 只 包含 抽象 的 坐标 
和 数据 。 为 了 直观 表现 ， 根 据 题目 所 给 数据 ， 通 过 二 维 三 次 多 项 式 插 值 的 
方法 ,使 用 MATLAB 工具 ， 大 致 绘制 出 八 种 重金 属 元 素 的 空间 分 布 图 和 区 
域内 的 海拔 等 值 线 图 ， 通 过 对 比分 析 重 金属 的 空间 分 布 。 关 于 每 个 区 的 污 
染 评 价 ， 以 《土壤 环境 质量 标准 》(GB15618-1995) 和 《绿色 食品 产地 环境 
质量 现状 评价 纲要 》( 试 行 ) (1994) 为 基本 标准 ， 分 别 计算 出 每 个 区 的 单 
因子 污染 指数 和 尼 梅 罗 污 染 综合 指数 ， 结 合 两 种 指数 对 各 个 区 域 的 污染 做 
ЖІ 

问题 二 、 根 据 数 据 作 出 的 重金 属 空间 分 布 图 和 海拔 等 值 线 图 ， 分 析 污 
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染 物 运 移 的 目 然 原因 和 人 为 原因 。 根 据 各 个 区 域内 的 污染 评价 ， 找 到 污染 
物 的 来 源 ， 得 到 重金 属 污染 的 原因 。 

问题 三 、 虽 然 有 重金 属 的 空间 分 布 图 ， 但 是 这 种 插值 结果 并 未 经 过 检 
验 ， 所 以 不 能 作为 查找 污染 源 的 最 终 手 段 ， 只 能 作为 参考 。 通 过 图 看 到 有 
坚 重 金属 分 布 广泛 ， 根 据 分 布 于 海拔 的 关系 ， 建 立 赋 权 值 的 重心 选 址 模型 ， 
采用 穷 举 法 ， 得 出 重金 属 污染 源 的 大 概 位 置 ， 再 通过 分 布 图 的 分 析 ， 进 行 
数据 筛选 ， 使 用 重心 坐标 选 址 模型 修正 污染 源 位 置 ， 有 些 重 金属 分 布 集中 ， 
直接 进行 数据 筛选 使 用 重心 坐标 选 址 模型 得 出 污染 源 位 置 。 

问题 上 四、 对比 三 次 多 项 式 插值 与 元 里 金 插 值 ， 对 比 污染 指数 评价 和 模 
糊 综合 评价 ， 对 比重 心 选 址 模型 与 反 演 模型 ， 通 过 对 比 ， 得 出 模型 的 优 缺 
点 。 如 果 搜 集 到 不 同 深度 的 土壤 检测 数据 ， 可 以 建立 迁移 模型 ， 得 出 污染 
的 时 间 ， 从 而 得 出 污染 是 由 污染 物 的 传播 还 是 由 于 污染 物 的 排放 引起 的 ， 
从 而 确定 污染 源 的 位 置 。 


3 模型 假设 
I 假设 该 城区 土壤 PH 值 为 6.5~7.5。 
П 不 考虑 该 城区 不 同 区 域 土质 的 差别 。 
II 在 但 找 污染 源 确切 位 置 时 ， 可 以 忽略 浓度 分 布 较 小 的 地 区 的 数据 。 
IV 假设 每 处 的 土壤 可 溶性 及 体积 相同 ,在 计算 重心 时 只 考虑 浓度 即 可 。 
У 在 用 赋 权 值 的 重心 坐标 选 址 法 的 过 程 中 ， 筛 选 数据 之 后 ， 由 于 有 海 
拔 作为 权 值 ， 所 以 可 以 不 用 考虑 横向 作用 力 对 污染 源 位 置 的 影响 。 

4 符号 说 明 


表格 1---- 符 号 说 明 


Рі =1233......8) 土壤 中 编号 为 i 的 重金 属 污 染 物 
的 污染 指数 

c=123....8) 土壤 中 编号 为 i 的 重金 属 污 染 物 
的 实测 浓度 

544-1224...) 编号 为 i 的 重金 属 污染 物 的 评价 
标准 〈 背 景 值 ) 

Т) 土壤 中 重金 属 污 染 物 的 综合 污 

染指 数 
区 域 的 横 坐标 
区 域 的 纵 坐 标 
5 模型 的 建立 与 求解 


5. 1 问题 一 : 重金 属 元 素 的 空间 分 布 及 污染 程度 分 析 
5. 1. 1 重金 属 元 素 的 空间 分 布 
插值 法 是 函数 逼近 的 一 种 重要 方法 , 三 次 多 项 式 插值 主要 讨论 了 给 


定 采 样 点 处 的 函数 值 后 ， 构 造 一 个 三 次 插值 多 项 式 的 方法 。 


建立 三 次 多 项 式 插值 模型 ， 使 用 MATLAB 软件 [1] 代码 见 附录 (1)) 


绘制 出 该 城区 不 同 重金 属 的 二 维 分 布 图 和 海拔 等 值 线 图 [2]。 
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1---- 海 拔 等 值 线 图 
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2----As 的 分 布 


As 的 分 布 范围 为 1.61~ 30.13hg/g， 空 间 分 布 不 连续 ， 高 值 区 域 呈 现 
定 的 扩散 趋势 ， 并 且 图 中 左 半 部 分 的 扩散 现象 明显 [3]。 对 比 等 值 线 图 ， 发 


现 扩散 区 域 地 势 较 低 很 可 能 是 由 于 污染 物 的 运 移 造 成 的 污染 。 
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3----Cd 的 分 布 图 


Cd 的 含量 范围 为 40. 00~ 1619. 80ng/g， 空 间 分 布 趋势 不 明显 ， 高 值 区 
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域 呈 现 一 定 的 扩散 趋势 。 对 比 等 值 线 图 ， 发 现 扩 散 区 域 地 势 较 低 很 可 能 是 


由 于 污染 物 的 运 移 造 成 的 污染 。 
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图 表 4----Cr 的 分 布 图 
Cr 的 含量 范围 为 15. 32~920. 84hg/g, 空间 分 布 特征 明显 ， 高 值 区 集中 


在 {(х, у) хе [2500,5000] ye [4000,600 叶 区域。 污染 物 没 有 明显 扩散 。 
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5----Си 的 分 布 
Cu 的 含量 范围 为 2. 29~2528. 48ug/g， 空 间 分 布 特征 明显 ， 高 值 区 集中 


在 xe[0,5000],ye[2000,4000] 区 域 ， 分 布 相对 集中 。 
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6----Hg 的 分 布 图 


Hg 的 含量 范围 为 8.57 ~16000.00 ng/g， 空 间 分 布 不 连续 ， 高 值 区 集 
中 在 
{(x, y)| x є [0,5000| y є |4000,6000)) 
{(x, у) x €[12500,15000], y e[1000,3000] 
{(x, y)| x €[12500,17000], у є [8000,10000]} 


区 域 ， 基 本 不 分 散 。Hg 进入 土壤 后 95% 以 上 能 迅速 被 土壤 固定 ， 因 此 Hg 容 


易 在 土壤 中 积累 。 
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7----Ni 的 分 布 图 


Ni 的 含量 范围 为 4. 27~142. 50hg/g ， 空 间 分 布 特征 不 明显 ， 高 值 区 集 


HE f(x, у) х є [2500,5000] ys[4000.600 趾 区域， 但 分 散 比 较 均 匀 。 
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8----Pb 的 分 布 图 


Pb 的 含量 范围 为 19. 68-472. 48ug/g ， 空 间 分 布 特征 明显 , 高 值 区 集中 
ТЕ(6,у)|хе1|0,500| у е|р000,6000))Я1 {(x, y)| x є [2500,7500], у є [4000,6000] DX 4K, 在 


周边 地 区 扩散 。 
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图 9----Zn 的 分 布 图 
Zn 的 含量 范围 为 32. 86~3760. 82ug/g ， 空 间 分 布 特征 明显 ， 高 值 区 集 
中 在 {(x,y)|xe[7500,10000| у є [3000,5000] 和 
(x, y)| x e[12500,15000], ye[8000,1000 叶 区域， 扩散 较 少 。 
5. 1. 2 重金 属 元 素 污染 程度 分 析 
考察 人 类 活动 对 污染 情况 的 影响 ， 对 该 城区 各 功能 区 各 种 重金 属 污染 
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现状 的 评价 可 以 运用 单 因子 污染 指数 法 和 尼 梅 罗 污 染 综 合 指数 法 结合 来 完 


成 [4,5]。 


单 因 子 污染 指数 公式 
P= (1) 


尼 梅 罗 污 染 综 合 指数 公式 


(9 平均 十 (一 ) “最 大 


Ш (2) 
2 


为 表述 方便 ， 将 各 重金 属 元 素 编号 。 


表格 2---- 各 重金 属 元 素 编 号 


1 1 2 3 4 5 6 8 


As Cd CT Cu Hg Ni Pb Zn 


根据 《土壤 环境 质量 标准 》(GB15618 一 一 1995) 和 《绿色 食品 产地 环 


境 质量 现状 评价 纲要 》( 试 行 ) (1994) 确定 土壤 重金 属 单 因子 污染 指数 评 


价 和 尼 梅 罗 污染 综合 指数 评价 的 污染 等 级 划分 。 


表格 3---- 土 壤 重 金属 污染 等 级 划分 标准 [4] 


等 级 综合 污染 指数 污染 等 级 污染 水 平 
1 Ps < 0.7 安全 清洁 
2 0.7 < № <1.0 警戒 级 尚 清洁 
3 1.0 < Ps <2.0 轻 度 污染 土壤 受到 轻 度 污染 
4 2.0 < Py, 53.0 中 度 污染 土壤 受到 中 度 污 染 
5 В > 3.0 乍 污染 土壤 污染 已 相当 严 
重 


利用 Microsoft Excel 软件 对 实验 数据 进行 处 理 ， 得 到 各 功能 区 重金 属 
单 因 子 污染 指数 P 和 尼 梅 罗 污 染 综合 指数 P,， 结 合 表 2 可 以 得 到 各 功能 区 
的 环境 污染 级 别 如 下 表 。 


表格 4---- 综 合 污染 评价 结果 表 


整个 城 公园 绿 
功能 区 | 区 | 生活 区 | 工业 区 | 山区 | 交通 区 | 地 区 
As 1.8367 | 17417 | 2.0139 | 11222 | 1.5861 1.7389 
污染 等 级 | 轻 度 污 | 轻 度 污 | 中 度 污 | 轻 度 污 | 经 度 污 | 轻 度 污 
Са | 2.6739 | 2.2305 | 3.0239 | 11717 | 2.7693 2.158 
污染 等 级 | 中 度 污 | 中 度 污 | 重 污染 | 轻 度 污 | 中 度 污 | шы 
CT 1.9797 | 22265 | 1.7229 | 1.2568 | 1.8726 | 1.4077 
污染 等 级 | 轻 度 污 | 中 度 污 | 轻 度 污 | 轻 度 污 | 经 度 污 | 轻 度 污 
Cu 7.4907 | 3.7424 | 9.6621 13121 | 4.7129 | 22871 
重 污染 | 重 污染 | 重 污染 | 轻 度 污 | 重 污染 | 中 度 污 
污染 等 级 染 染 
Hg 14.0589 | 2.6583 | 18.3531 | 1.1703 | 12.7663 | 3.2854 
重 污染 | 中 度 污 | 重 污染 | 轻 度 污 | 重 污染 | 重 污染 

污染 等 级 染 染 
Ni 1.5121 | 1.4911 1.6106 | 1.2561 1.4325 1.2431 


轻 度 污 | 轻 度 污 | HEY | HE | 轻 度 污 | 轻 度 污 


污染 等 级 | 染 | 染 | | 染 | 染 ра 
Pb 2.5835 | 2.2294 | 3.0013 1.1794 | 2.0494 1.9584 

中 度 污 | 中 度 污 | 重 污染 | 轻 度 污 | 中 度 污 | 轻 度 污 
污染 等 级 | Е жж 
Zn 3.4915 | 3.4349 4.028 1.0622 | 3.5196 | 2.2354 


重 污 染 | 重 污 染 | 重 污染 | HEH | 重 污染 | 中 度 污 


污染 等 级 染 染 
整个 城区 3.1704 | 13.5331 | 1.2532 | 94264 | 2.7343 
综合 指数 


EIR | 重 污染 | HEY | 重 污染 | 中 度 污 


污染 等 级 染 染 


由 表 3 可 知 : 

就 整个 城区 而 言 As、Cr、Ni 为 轻 度 污染 ，Cd 和 Pb 为 中 度 污染 Си. 
Hg 和 Zn 属于 重 污 染 ， 其 中 Cu 和 Не 的 污染 尤为 严重 ,: 

生活 区 As 和 Ni 为 轻 度 污染 , Са. Сг. Не 和 Pb 为 中 度 污 染 ，Cu 和 Zn 


属于 重 污染 ， 整 个 功能 区 评定 为 重 污染 ; 


工业 区 Сі. Ni 为 轻 度 污染 ，As 为 中 度 污 染 ，Cd、Cu、Hg、Pb、Zn 
为 重 污 染 ， 其 中 Cu 和 Hg 的 污染 特别 严重 ， 整 个 功能 区 评定 为 重 污染 ; 


山区 全 部 为 轻 度 污染 ， 土 壤 状 况 相 对 恨 好 ， 整 个 功能 区 评定 为 轻 度 污 


交通 区 Ав. Ог. Ni 为 轻 度 污染 ，Cd 和 Pb 为 中 度 污 染 ，Cu、Hg 和 Zn 
为 重 污 染 ， 其 中 Hg 的 污染 最 为 严重 ， 整 个 功能 区 评定 为 重 污 染 ; 

公园 绿地 区 的 As, Gr, Ni, Pb 为 轻 度 污染 ，Cd、Cu、Zn 为 中 度 污 染 ， 
Hg 为 重 污染 ， 整 个 功能 区 评定 为 中 度 污 染 。 

所 以 从 综合 评价 的 结果 来 看 ， 该 城区 的 重金 属 污染 较为 严重 ， 尤 其 是 
工业 区 和 交通 区 ， 情 况 不 容 乐观 。 


5. 2 问题 二 : 分 析 重 金属 污染 的 主要 原因 

通过 分 析 对 比重 金属 分 布 图 和 海拔 等 值 线 图 可 以 猜测 各 重金 属 元 素 大 
致 的 来 源 。 

As、Cd 既 有 集中 分 布 的 区 域 ， 又 有 广泛 扩散 的 区 域 ， 而 且 广 泛 分 散 分 
布 在 海拔 较 低 的 部 分 ， 根 据 推 流 运 移 理论 ， 猜 测 这 一 部 分 是 由 空气 或 者 流 
水 传播 过 来 的 。 故 在 高 处 浓度 集中 分 布 的 地 区 ， 可 近似 划 定 为 污染 的 源头 
所 在 区 域 。 

Cr, Cu, Hg, Ni 只 有 集中 分 布 的 区 域 ， 说 明 仪 在 这 一 个 地 带 富 集 ， 说 
明 这 一 区 域 为 富 集 区 或 源头 所 在 地 。 

Pb, Zn 虽然 既 有 集中 分 布 的 区 域 ， 又 有 广泛 扩散 的 区 域 ， 但 是 与 海拔 
关系 不 大 ， 所 以 扩散 视 为 人 为 因素 。 

再 根据 污染 评价 的 结果 对 每 个 区 作 分析 。 因 评价 时 选用 了 该 城区 上 自然 区 
的 重金 属 元 素 含量 为 背景 值 ， 所 以 不 必 考 虑 目 然 原 因 对 其 来 源 的 影响 。 

生活 区 As 和 Ni 为 轻 度 污染 ,Cd、Cr、Hg 和 Pb 为 中 度 污 染 ，Cu 和 Zn 


属于 重 污 染 ; 


工业 区 Cr, Ni 为 轻 度 污染 ，As J HPs, Cd, Cu, Hg, Pb, Zn 
为 重 污染 ， 其 中 Cu 和 Heg Mis HR EE 

山区 全 部 为 轻 度 污染 ， 土 壤 状 况 相 对 良好 ; 

交通 区 As、Gr、Ni 为 轻 度 污染 ，Cd 和 Pb 为 中 度 污 染 ，Cu、Hg 和 Zn 
为 重 污染 ， 其 中 Hg 的 污染 最 为 严重 ; 

公园 绿地 区 的 As, Gr, Ni, Pb 为 轻 度 污染 ，Cd、Cu、Zn 为 中 度 污 染 ， 
Hg 为 重 污染 。 

由 所 给 这 些 数据 , 分 别 从 重金 属 空间 分 布 图 和 污染 评价 两 个 角度 综合 分 
析 并 结合 重金 属 用 途 [6] 得 到 如 下 结论 。 


表格 5---- 各 功能 区 重金 属 污染 主要 原因 


化 妆 品 ` 照 明灯 、 等 含有 大 量 的 Нр, 
燃 煤 排 放大 量 的 Hg 和 Zn; Ж 
电池 中 Cd 量 较 高 ; 劣质 化 妆 品 中 含 
有 大 量 的 Сг; Pb 主要 来 源 于 各 种 油 
生活 区 BW Mh Ei, IW, hs, 
DIRS Н/Ж. Лат. YAH Ж 


PWR, ER, WEEE WEE hs 
自来水 管 等 。 这 些 重金 属 元 素 有 些 
在 垃圾 堆放 过 程 中 经 过 自然 目 然 沉 
降 和 两 淋 沉 降 进 入 土壤 ， 有 些 通过 
食物 链 进入 人 体 ， 排 汇 物 含有 大 量 
的 重金 属 元 素 ; 有 些 随 粉 全 进 入 十 


BE 


工业 区 


Cu EBERLE MED 的 开采 和 论 
炼 、 金 属 加 工 、 机 械 制 造 、 钢 铁 生 
产 等 ，Hg 主要 来 源 于 仪表 三 、 食 盐 
HE, SHENK, РЬ 的 主要 来 
源 为 电镀 ; Cd 要 来 源 有 电镀 、 采 矿 、 
冶炼 、 燃 料 、 电 池 和 化 学 工业 每 ; 
Zn 的 主要 来 源 开采、 冶炼 加 工 、 机 
械 制 造 以 及 镀 笑 、 仪 器 仪表 、 有 机 
会 合成 和 造纸 等 工业 的 排放 ; As 主 


要 来 源 于 采矿 、 治 金 、 化 学 制药 等 


РЬ. Zn, Cd, Cu 主要 来 源 于 含 Pb 
汽油 的 燃烧 和 汽车 轮胎 磨损 产生 
的 粉尘 ， 经 过 自然 沉降 和 雨 淋 沉 


降 进 入 土壤 


农药 、 化 肥 和 地 膜 是 重要 的 农 


用 物资 ， 但 长 期 不 合理 施用 可 能 4 
致 土壤 重金 属 污 染 ， 个 别 农药 中 含 
有 Нр. Аѕ. Си. 20 等 重金 属 。 大 
部 分 Hg 随 消毒 种 子 进入 土壤 ， 含 
公园 绿地 区 As 农药 的 施用 造成 As 在 土壤 中 的 
积累 ， 杀 菌 农 药 造成 Cu、Zn 累积 ， 
随 磷肥 进入 土壤 的 Cd 也 一 直 受 到 
关注 。 而 且 污水 灌溉 使 公园 绿地 区 


受 工 业 区 重金 属 含量 的 影响 ，Cd、 


Cu, Не. Zn 的 也 逐渐 积累 ，Pb 的 
迁移 性 弱 , 所 以 灌溉 区 Pb 的 累积 特 


点 是 离 污染 源 近 土 壤 中 含量 高 。 


5. 3 问题 三 ， 分析 重 金属 污染 物 的 传播 特征 ， 确 定 污染 源 的 位 置 
5. 3. 1 金属 污染 物 的 传播 特征 

对 比重 金属 分 布 图 和 海拔 等 值 线 图 可 以 得 到 : 

Ass Cd 既 有 集中 分 布 的 区 域 ， 又 有 广泛 扩散 的 区 域 ， 并 且 广 泛 分 散 分 
布 在 海拔 较 低 的 部 分 ， 根 据 推 流 运 移 理论 ， 猜 测 这 分 散 分 布 的 部 分 是 由 空 
气 或 者 流水 传播 过 来 的 。 

Cr, Cu, Hg, Ni 只 有 集中 分 布 的 区 域 ， 说 明 这 些 重金 属 仅 在 这 一 个 地 
带 富 集 ， 说 明 这 一 区 域 为 富 集 区 或 源头 所 在 地 。 

РЬ. Zn 虽然 既 有 集中 分 布 的 区 域 ， 又 有 广泛 扩散 的 区 域 ， 但 是 与 海拔 


关系 不 大 ， 所 以 扩散 视 为 人 为 因素 。 
5. 3. 2 确定 污染 源 的 位 置 
先 建立 带 权 值 的 重心 坐标 选 址 法 [7] 确 定 As、Cd 污染 源 的 大 概 位 置 。 
带 权 值 的 重心 坐标 选 址 法 , 是 指 基于 污染 物 海拔 高 度 的 不 同 附 上 不 同 的 
权 值 ， 然 后 用 重心 坐标 选 址 的 方法 找 出 污染 源 的 位 置 。 由 于 重心 坐标 的 位 
置 是 所 有 重金 属 元 素 相对 于 其 位 置 重力 抢 最 小 的 位 置 , 用 MATLAB 采用 穷 
举 法 的 算法 进行 运算 ， 找 到 重心 距 之 和 最 小 的 位 置 。 设 重心 坐标 的 位 置 为 
(xz，y%)， 原 题 所 给 数据 中 任意 一 处 的 元 素 坐 标 为 《”，yw，)， 浓 度 为 ci， 
体积 为 V,。 设 每 个 采样 点 的 重力 矩 为 F(x,y)。 


F(x, y) = V, хс, хуб, —X)° +09, 0) (3) 
Обу) = 3) F(x,.9,) (4) 


然后 利用 穷 举 法 使 用 MATLAB THRE AUG LE 2), (322. 
解 得 : As 的 重心 坐标 为 〈16800,10080)，Cd 的 重心 坐标 为 (15600,9520)。 
结合 图 形 发 现 所 求 的 As 的 重心 坐标 接近 富 集 区 为 
(о, ух е115140,19767), у €[8287,11933]}, 
Cd 的 重心 坐标 接近 富 集 区 为 
{(x, y)}x є[2486,4948], у €[4603,7293]} , 
观察 其 他 重金 属 元 素 的 分 布 ， 并 无 出 现 扩散 ， 但 是 出 现 多 位 置 的 富 集 ， 
分 别 筛选 数据 ， 选 取 富 集 区 以 内 的 采样 点 。 
表格 6---- 各 种 重金 属 富 集 区 域 


人 人、 
重 HE x min * max У тіп y max 


属 
Ав 15140 19767 8287 11933 
Са 2486 4948 4603 7293 
Ст 2586 4777 4357 6018 
Са 2383 3927 2110 4357 
Hg 2383 3927 2110 3971 
Ni 3299 4948 4357 7293 
Pb 3526 6924 4080 5696 
Zn 13797 14325 8666 11941 


再 次 进行 重心 坐标 选 址 的 MATLAB 人 穷 举 算法 《代码 见 附录 C2), 032), 
BSR, hl 


作 如 下 表格 。 
表格 7---- 污 染 源 位 置 坐标 
重金 属 元 素 х(т) y(m) 
As 18100 10080 
Cd 3900 5950 
Cr 3800 5320 
Са 2400 3710 
Hg 2700 2310 
Ni 3300 6020 
Pb 4800 4830 


Zn 13800 9590 


5. 4 问题 四 : 模型 评价 、 改 进 及 新 模型 引入 


5. 4. 1 模型 的 优点 

本 文 分 析 重 金属 空间 分 布 时 运用 了 二 维 三 次 多 项 式 插值 的 方法 绘制 重 
金属 分 布 图 ， 多 项 式 插 值 是 根据 已 知 采 样 点 的 数据 值 拟 合 一 个 多 项 式 函 数 ， 
用 这 个 多 项 式 函 数 来 表达 连续 分 布 的 地 理 空间 的 方法 ， 这 种 方法 简单 、 使 
用 方便 。 

运用 单元 型 污染 指数 和 尼 梅 风土 壤 污 染 综合 指数 评价 法 结合 对 该 城区 
污染 情况 做 出 了 分 析 ， 单 元 型 污染 指数 法 单独 对 各 个 功能 区 每 种 金属 的 污 
染 进 行 分 析 ， 内 梅 罗 土壤 污染 综合 指数 法 对 整个 城区 每 种 元 素 的 综合 污染 
进行 分 析 ， 数 学 过 程 简 捷 ， 物 理 概 念 清晰 ， 便 于 决策 者 做 出 综合 决策 。 

在 问题 二 中 分 析 污 染 来 源 时 结合 迁移 理论 和 分 布 图 ， 结 果 更 准确 。 

问题 三 确定 污染 源 位 置 时 ， 先 进行 数据 筛选 再 运用 重心 选 址 模型 ， 避 免 
了 大 量 的 运算 ， 并 且 此 种 方法 要 求 各 样本 点 固定 ， 适 用 于 本 题 。 


5. 4. 2 模型 的 缺点 

尽管 运用 多 项 式 插 值 计算 简 单 ， 但 这 种 算法 误差 较 大 。 

尼 梅 罗 土 壤 污 染 综 合 指数 评价 法 过 于 突出 最 大 污染 因子 对 土壤 污染 的 
影响 ， 没 有 考虑 到 权重 因素 ， 没 有 把 危害 性 最 大 的 因子 放 在 重要 角度 进行 
分 析 。 
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确定 污染 源 的 位 置 除 运用 了 重心 选 址 模型 , 在 函数 图 形 不 均匀 的 情况 下 
会 出 现 结果 突变 现象 ， 影 响 污染 源 的 确定 。 


5. 4. 3 模型 的 改进 


(0 “对 于 多 项 式 插值 算法 误差 大 的 问题 ， 我 们 可 以 运用 克 里 金 插值 进 
行 改进 。 

克 里 金 插值 法 假设 区 域 空间 变量 x 在 空间 上 的 分 布 满足 函数 f(x), 并 将 
x 的 空间 分 布 分 解 为 三 部 分 : 变量 x 在 空间 上 分 布 的 主体 趋势 ; 变量 x 具有 
一 定 相 关 性 的 随机 影响 变量 ， 这 些 随 机 变量 受 x 分 布 的 影响 ， 随 机 噪声 项 
与 x 的 空间 分 布 不 具备 相关 性 ， 且 与 x 分布 的 整体 趋势 无 关 。 

OO: (5) 
9, р(х) Ира x 分 布 的 整体 趋势 函数 ，wlx) 是 随机 影响 变量 ，< 是 
随机 噪声 项 。 

克 里 金 插值 法 又 称 为 空间 自 协 方差 最 佳 内 插 法 ， 不 仅 可 以 获得 预测 结 
果 ， 还 能 够 获得 预测 误差 ， 有 利于 评估 插值 结果 的 不 确定 性 。 克 里 金 插值 
法 是 精度 较 高 的 一 种 算法 ， 但 这 种 算法 运算 量 大 ， 计 算 耗费 时 间 比 较 多 ， 
这 里 我 们 不 做 详细 计算 。 


2) 为 了 更 合理 地 反映 土壤 环境 的 污染 程度 ， 应 考虑 危害 性 最 大 的 重 
金属 的 权重 因素 。 改 进 的 梅 罗 综合 污染 指数 计算 公式 为 : 
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0, | (6) 
其 中 ,为 权重 值 最 大 的 污染 因子 6 与 s, 的 比值 。 
各 种 重金 属 对 环境 和 人 体 的 危害 程度 有 很 大 差别 ， 一 般 的 ， 某 种 重金 
属 的 危害 性 与 
该 污染 因子 的 背景 值 基本 上 呈 反 比例 的 内 在 关系 ,首先 将 各 重金 属 s, 的 背景 
值 由 小 到 大 的 顺序 排列 ， 并 将 其 最 大 值 ,和 ,比较 ， 令 R, 表 示 第 i 中 重金 


属 的 相关 性 比值 。 


R,= зы (7) 
Aw, ЕЕ 8 (8) 
Ук 


i=l 


其 中 ，w 为 第 i 种 污染 因子 的 权重 值 ， 显 然 》w =1。 


通过 简单 计算 可 以 得 到 各 重金 属 的 权重 值 如 下 表 。 


表格 8---- 各 重金 属 权重 值 


Was Уса We, Woy Wy 


6 Wi Wpp Му, 


0.0076 | 0.2092 | 0.0009 | 0.0020 | 0.7768 | 0.0022 | 0.0009 | 0.0004 


HR ”可 以 看 出 权重 值 最 大 的 重金 属 是 Hg， 污 染 等 级 的 划分 同 表 2, 
根据 采样 点 重金 属 的 浓度 ， 运 用 公式 〈6) 计算 出 各 功能 区 的 尼 罗 梅 综合 污 
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染指 数 如 下 表 。 


表格 9---- 改 进 后 的 尼 梅 罗 综 合 污 染 评 价 指数 


生活 
功能 区 区 | 工业 区 山区 | 交通 区 公园 绿地 区 
综合 污染 指 | 3.2004 | 18.3531 | 1.2412 | 12.7663 3.2854 
数 
重 污 | 重 污染 | HER | Risk 重 污 染 
污染 级 数 ik 


对 比 表格 3 和 表格 8 发 现 各 功能 区 的 综合 污染 指数 有 所 提高 , 工业 区 和 
交通 区 的 提高 尤其 大 ， 这 是 因为 改进 的 尼 梅 罗 评 价 法 考虑 了 危害 性 最 大 的 
重金 属 Hg 的 权重 因素 ， 而 由 Hg 在 工业 区 和 交通 区 的 污染 来 源 最 广泛 可 知 


改进 后 的 评价 方法 更 科学 合理 。 


(3) 


为 了 克服 重心 选 址 模型 出 现 的 突变 现象 ， 我 们 可 以 运用 反 演 模型 
进行 改进 。 反 演 是 把 空间 所 有 力 一 向 同时 反 过 来 ， 利 用 各 样本 点 的 重金 属 
浓度 ， 对 污染 源 的 位 置 和 强度 进行 反 演 ， 可 以 得 到 概率 分 布 情况 ， 对 反 演 


后 的 统计 量 进行 对 比 ， 进 而 分 析出 污染 源 的 位 置 。 


5. 4. 4 收集 信息 引进 模型 


为 研究 城市 地 质 环境 的 演变 模式 ， 还 应 收集 信息 解决 对 应 问题 。 
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() ЖЛ ЕЛ, MRS VTE Ж. Жа. ПАУ ЛУ, 
FE LEEPER TE ЛАЛА, EB AER RSE ЕЕ ЖЕЛЕ ЖИА ҒА АЗЫ 
不 同 而 不 同 。 其 中 PH 值 是 决定 土壤 重金 属 存在 形态 的 关键 因素 ， 因 此 收集 
不 同 采 样 点 的 PH 值 ， 根 据 国家 环境 质量 评定 标准 ， 用 单 因子 污染 指数 和 尼 
梅 罗 土壤 污染 综合 指数 法 计算 出 各 个 金属 元 素 更 精确 的 污染 级 别 。 


2 由 于 大 多 数 重 金属 污染 物 在 土壤 中 的 衰减 速度 和 转化 过 程 比 
较 慢 ， 污 染 物 的 总 质量 和 总 浓度 变化 可 以 忽略 不 计 ， 随 着 时 间 的 推移 ， 
己 进 入 土壤 中 的 重金 属 污染 物 在 气流 或 者 水 流 作 用 下 产生 空间 位 置 上 
的 转移 ， 致 使 不 同 土 层 中 污染 物 的 浓度 发 生变 化 。 推 流 迁移 只 改变 污染 
物 所 处 位 置 ， 并 不 改变 污染 物 的 质量 和 浓度 。 因 此 利用 数学 抽象 迁移 通 
量 (单位 时 间 通 过 单位 面积 的 量 ) 来 描述 污染 物 的 推 流 迁移 。x、y 和 z 
方 癌 的 污染 物质 量 通 量 为 


Ат, -и,с (9) 
Ат, -и,с (10) 
Am. = ис (11) 


式 中 : Am、Am、Am.---- 在 x、y 和 z 方 向 上 污染 物 推 流 迁 移 通 量 ; 
их uu,、u.---- 在 XxX、y 和 z 方 向 上 环境 介质 流速 分 量 ; 
c---- 污 染 物 在 环境 介质 中 的 浓度 。 
因此 为 研究 城市 地 质 环 境 的 演变 模式 , 可 收集 每 个 采样 点 对 应 的 深度 、 


土壤 质地 背景 、PH 值 和 土壤 含水 率 [8]， 运 用 推 流 迁 移 模 型 可 以 计算 出 不 
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同样 本 点 的 污染 物 浓度 ， 进 一 步 可 以 分 析出 污染 时 间 和 污染 源 的 位 置 。 
6 参考 文献 

[1] 周 建 兴 AXP BAER ARE, (MATLAB 从 入 门 到 精通 》, 北京 : 人 

民 邮 电 出 版 社 ，2008. 11 (2010. 1 重印 )。 

[213852 XA, (615 与 空间 分 析 一 一 原理 与 方法 》， 北 京 : 科学 出 版 社 ， 

2006 

[3] 许 金条 , 《空间 分 布 特征 与 相关 性 》, 《基于 GIS 的 金华 市 区 土壤 中 金属 

含量 空间 分 布 与 污染 评价 研究 》，20080501，35-37, 2008 年 

[4] ЕЕ, 《环境 评价 》 上 海 : 同济 大 学 出 版 社 ，1999. 9。 

[5] 肖 再 亮 ， 兰 石 ， 田 犀 ， 陈 钢 ， 苏 培 生 , 《 利 州 区 农产品 基地 土壤 重金 属 


污染 评价 》， 环 境 评 价 ， 第 29 卷 第 3 期 : 131-135，2010 年 。 

[6] 刘 青松 , 《环境 污染 与 防治 技术 》 北京 : 中 国 环境 科学 出 版 社 ，2004 年 
3 月 。 

[7] 苗 兴 东 ” 李 映 红 ” 范 存 军 ,《 中 心 选 址 探讨 ),《 交 通 标准 化 》， 第 134 期 : 
50-52, 2004 年 。 

[8] 王 洋 ХІ ЖО, CELE PH 值 对 冻 融 黑土 重金 属 锌 赋 存 形态 的 影 
响 》， 中 国 科学 院 东 北 地 理 与 农业 生态 研究 所 ， 长 春 ，130012。 


25 


附录 
(1) 元 素 分 布 代码 
>>х-0(,2:2); 
>>у-0(:,3:3); 
>> 7=0)(:,6:6)'; 
>> х1=0:100:30000; 
>> у1=0:70:21000; 
>> [Х,Ү|=теѕһогпа(х1,у1); 
>> H=griddata(x,y,z,X, Y,'cubic'); 
>> pcolor(X, Y,H);shading flat 
>> colorbar; 
xlabel(' 横 坐标 х(т)) 
ylabel(' 纵 坐标 у(т)) 
title('As 的 空间 分 布 ') 
(2) 单个 功能 区 选 址 代码 
х-Ав(:,2:2); 
у-А5(:,3:3); 
72=А5(:,6:6)'; 
n=length(x); 
іх=х; 
їу=у; 


(7-7; 
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min=inf; 

for bx=(0:100:30000); 

for by=(0:70:21000); 
[bx,by |] =meshgrid(bx,by); 
?-0; 

for j=1:1:n; 

M=t2(J); 
T=(sqrt((bx-tx(j))."2+(by-ty(j)).M2)); 
?-27--М.%Т; 

епа 

if min>Z 

min=Z; 

x=bx; 

у=Ыу; 

епа 


епа 


(3) 整个 城区 选 址 代码 
x=0(:,2:2)'; 
у=0(:,3:3)"; 
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7=0(:,6:6)'; 
h=Q(:,4:4)'; 
n=length(x); 

іх=х; 

їу=у; 

(7-7; 

th=h; 

min=inf; 

for bx=(0:200:30000); 
for by=(0:140:21000); 
[bx,by ]=meshgrid(bx,by); 
?-0; 

for j=1:1:n; 

M=t2(J); 
T=(sqrt((bx-tx()).42+(by-ty(j)).*2)); 
P=th(j); 

Z=Z+M.* T.*P; 

end 

if min>Z 

min=Z; 

x=bx; 


y=by; 
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(4) Excel 示例 元 素 As I) fifi 26 0) 


编号 X Y 海拔 | 区 域 | 浓度 
30 17904 | 8287 29 1 5.83 
29 16301 | 8299 24 1 5.41 
237 19072 | 8519 36 4 6.09 
21 15467 | 8658 17 1 2:91 
257 19767 | 8810 46 4 6.98 
25 16428 | 9069 20 1 4.38 
196 17144 | 9081 20 4 2:91 
186 15248 | 9106 16 4 2527 
212 18393 | 9183 26 4 1.96 
239 17198 | 9810 37 4 3.89 
245 18134 | 10046 41 4 30.13 
38 16289 | 10072 43 1 5.41 
222 15198 | 10100 28 4 7.34 
117 17044 | 10691 93 3 2.72 
274 17814 | 10707 64 4 6.35 
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96 18738 | 10921 53 10.99 
271 16267 | 11058 60 4.38 
72 15140 | 11101 30 3.37 
113 19007 | 11488 84 6.14 
115 18413 | 11721 88 3.69 
248 17087 | 11933 43 8.50 
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